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Abstract 

The construction of a computer code to calculate the cross sections for the 
spin-polarized processes e~7 — > e~7, e~77, e~~e + e~ to order-o 3 is described. 
The code calculates cross sections for circularly-polarized initial-state photons 
and arbitrarily polarized initial-state electrons. The application of the code to 
the SLD Compton polarimeter indicates that the order-a 3 corrections produce 
a fractional shift in the SLC polarization scale of —0.1% which is too small and 
of the wrong sign to account for the discrepancy in the Z-pole asymmetries 
measured by the SLD Collaboration and the LEP Collaborations. 
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I. INTRODUCTION 



For some years, polarized Compton scattering, the scattering of circularly-polarized pho- 
tons by spin-polarized electrons, has been used to measure the degree of polarization of one 
particle or the other. Circularly-polarized gamma-ray photons from nuclear decays have 
been polarization-analyzed by measuring the asymmetry in the rates of backscattering from 
magnetized iron foils (as the foil magnetization is reversed). Similarly, the polarization of 
electrons in high-energy storage rings and accelerators has been determined from scatter- 
ing asymmetries of the accelerated beam with beams of optical laser photons. Until quite 
recently, all such measurements have made use of tree-level (order-a 2 ) expressions for the 
polarized Compton scattering cross section. 

Part of the reason for this has been the unavailability of a next-to-leading-order calcula- 
tion that is packaged in an easily usable form. The first calculation of the order-a 3 virtual 
and real-soft-photon corrections to unpolarized Compton scattering was published by Brown 
and Feynman in 1952 This calculation was not confirmed until 1972 when Tsai, De- 
Raad, and Milton (TDM) published the same corrections for the polarized case 0. The 
TDM calculation, by itself, is sufficient to interpret the results of measurements involving 
longitudinally polarized electrons for which the presence of additional energetic photons in 
the final state can be excluded. This is often the case for measurements of gamma-rays that 
have been scattered from magnetized iron targets. However, accelerator-based polarimeters 
are often designed to measure transverse electron polarization and generally cannot dis- 
tinguish between single-photon and multiple-photon final states. These shortcomings were 
addressed in 1987 by Gongora and Stuart (GS) who published the matrix elements for the 
hard-photon corrections in a spinor-product form that is suitable for numerical evaluation 
0. Their publication also includes spinor-product expressions for the matrix elements of 
six gauge- invariant tensors used by TDM to calculate their result. These expressions permit 
the application of the TDM virtual corrections to the case of general initial- and final-state 
electron spin directions. Finally, in 1989, the complete set of virtual, soft-photon, and hard- 
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photon radiative corrections to polarized Compton scattering was calculated independently 
by Veltman Veltman's paper describes her calculation qualitatively and presents a re- 
sult in numerical form for two specific cases of an accelerator-based longitudinal polarimeter. 
Unfortunately, it does not include detailed expressions for the final result nor is the result 
checked against the TDM or GS calculations. 

One of the specific cases discussed by Veltman, the case of a 50 GeV longitudinally- 
polarized electron colliding with a 2.34 eV photon, is quite close to that of the SLD Compton 
polarimeter (a 45.65 GeV longitudinally-polarized electron colliding with a 2.33 eV photon). 
This polarimeter is a key component in the measurement of the left-right Z-boson production 
asymmetry A° LR which has been performed over several years by the SLD Collaboration 
0. At the current time, the measured value of A° LR is approximately 8% larger than 
the value of the comparable quantity extracted from measurements of six different Z-pole 
asymmetries by the four LEP Collaborations 0. Since the SLD and LEP measurements 
differ by approximately three standard deviations, the discrepancy is more likely to be 
due to systematic effects than to statistical fluctuations. One possible systematic effect is 
the absence of radiative corrections from the interpretation of the SLD polarimeter data. 
Veltman's calculation implies that radiative corrections would shift the SLD polarization 
measurements by —0.1% of themselves which is far too small and of the wrong sign to 
account for the 8% discrepancy. 

This paper describes a complete order-a 3 calculation of polarized Compton Scattering. It 
was undertaken primarily to check the calculation of Veltman and to determine if radiative 
corrections to Compton scattering could be responsible for discrepancy between the LEP 
and SLD measurements of Z-pole asymmetries. A second goal was to develop a computer 
code which could applied to a variety of present and future experimental situations. The 
main ingredients of this code, the TDM and GS calculations, are sufficient for all present day 
experimental situations. However, it is likely that a very high energy linear electron-positron 
collider will be constructed somewhere in the world in the coming decade. Polarized beams 
are planned for all of the designs now under discussion. All of these projects incorporate 

3 



Compton Scattering polarimeters into the optical designs of their final focusing systems. If 
these polarimeters use optical lasers, the e~7 center-of-mass energies will be above threshold 
for the production of final state e + e~ pairs. Since the process e~7 — > e~e + e~ occurs at order- 
ed, it must also be included in the computer code. A calculation of the matrix element for 
this process, based upon the techniques of Ref. [||], is described in Section [11 D[ 

The following sections of this paper describe the construction and operation of the 
Fortran-code COMRAD which calculates the order-a 3 cross section for polarized Comp- 
ton Scattering. Section || describes the ingredients of the calculations which the code is 



based. Section |TJ describes the actual implementation of the various calculations and sev- 
eral cross checks that were performed. Section [TV] describes the application of the code to 
several cases of interest. And finally, Section [V] summarizes the preceding sections. 



II. INGREDIENTS 

This section describes the ingredients used to construct the code COMRAD. The hard 
photon corrections, virtual photon corrections, soft photon corrections, and e~e + e~ cross 
sections are discussed in the following sections. Since the e~e + e~ cross section calculation 
makes use of the techniques used to calculate the hard photon corrections, some technical 



details are presented in Section |II A| that facilitate the description of the original work 
presented in Section [II D . 



A. Hard-Photon Corrections 

The calculation of the cross section for the process e~7 — > e~77 is based upon the matrix 
element calculation of Gongora and Stuart ||. Their calculation is the first application of 
numerical spinor product techniques |7[] to a case involving massive spinors. These techniques 
allow one to express any amplitude as a function of the scalar products of two massless 
spinors u±(p) and their conjugates u±(p). The subscripts refer to positive and negative 
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helicity states of a massless fermion of momentum p. The only two non-vanishing scalar 
products, 

S+(Pl,P2) = U + (p 1 )u^(p 2 ) = -S+(p2,Pl) (1) 

s-(pi,p 2 ) = U-{px)u + {p 2 ) = -s+(p 1 ,p 2 )*, (2) 

are easy to evaluate numerically. Gongora and Stuart define the photon polarization vector 
in terms of these quantities so that it is free of axial-vector components and can be used 
with massive currents, 

4(<z, q) = ± fK V ^ (g)V^(g). ( 3 ) 

where: the ± subscript refers to the helicity of initial-state photons (final-state photons 
have opposite helicities), q is the photon momentum, and q is an arbitrary massless vector. 
Massive spinors of arbitrary spin direction are defined in terms of massless spinors as follows, 

u(p, s) = S+iP ^ P2) u + ( Pl ) + U_(P2) (4) 
u(p, s) = - S ~^ 1,P2 ^ + (pi) + w-(p 2 ) (5) 

where m is the electron mass and the massless vectors, p\ and p%, are defined in terms of 
the momentum and spin vectors, p and s, as follows, 

Pi = ^ (P + ms ) ( 6 ) 
P2 = ^ (P ~ ms ) • ( 7 ) 

The actual calculation involves the evaluation of a single Feynman amplitude 
D\\'\"(q,q;q',q';q",q") for the process e~(s) — ► e _ (s')7(A)7(A')7(A") (shown in Fig. |I]) 
where: A, A', and A" label the helicities of the three photons (+ or — ); and s and s' are 
the spin-vectors of the initial- and final-state electrons, respectively. The matrix element 
■M-\;\'\"(s, s') for the process e~(s)7(A) — > e~(s / )7(A / )7(A") can then be constructed from 
the function -Daa'A" by reversing the momenta of single photons and by interchanging the 
momenta and helicities of the remaining identical photons, 



M X ;X'X"(s, s') = D X \>\»(-q, q; <?', q'\ q", q") + D X \»\>(-q, q; q", q"\ q\ <f) 
+ Dyx\"(q', q; -q, q; q", q") + Dyxx'(q", q"] -q, q\ q', q') 
+ Dyx'x(q", q"] q', q'\ -q, q) + Dyx"x(q', q'\ q" , q"\ -q, q), (8) 

where q, q', and q" are the momenta of the incident and final-state photons, respectively. 
Note that each of the six terms in Eq. |8] corresponds to an ordinary Feynman diagram. 

Two technical issues are relevant to the present discussion and to the presentation of 
Section |11 D|. The first issue concerns the choice of the arbitrary massless momenta: q, 



q', and q" . Gongora and Stuart point out that one can substantially simplify some of the 
expressions by a judicious choice of these auxiliary momenta. They present results for two 
equivalently-simple sets of momenta. This approach provides an important cross check (one 
must find identical results for both sets of auxiliary momenta) and greatly facilitated the 
debugging of the GS manuscript and the computer code. A number of typographical errors 
were discovered in Ref. and are listed in Appendix [A|. One should note that the first set 
of auxiliary momenta (used to calculate GS Eqs. 3.3.1-3.10.2) always produces singularities 
when the initial state electron is longitudinally polarized whereas the second set (used to 
calculate GS Eqs. C.l.l-C.8.2) never develops singularities so long as both photons have 
non-zero energy. 

The second issue concerns the evaluation of spinors of negative momenta. In order to 
preserve the following (very useful) relationship, 

-(l±7s) 3^=«±(p)«±(p), (9) 

it is necessary to define negative momentum spinors in the following manner, 

u±(-p) = iu±(p), u±(-p) = iu±(p). (10) 

This, in turn, implies that spinor products of negative arguments behave as follows, 

s±(~qi, q*) = s±(q u -q 2 ) = is±(q u q 2 ) (11) 
s±(-9i, -92) = s±(qi, q 2 ), (12) 
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and that external photon polarization vectors are invariant under the transformation q — > —q 
(see Eq. |), 

e"(-g,<7) = e"( g ,g). (13) 

The actual cross section for the process e~(s) r y(X) — > e~77 is calculated in the center-of- 

mass (cm) frame from the matrix element given in Eq. || using the following expression, 

d 5 o- (1) 1 2 

dW^ 7 (a ' A) = U(2^E,(E e + P e ) TfJMwvMl , (14) 

where: _E e and P e are the energy and 3-momentum of the incident electron, E 1 is the energy 
of the incident photon, E' e and Q' e are the energy and direction of the final state electron, 
E' is then energy of one of the final state photons, and 0^ is the azimuth of the final state 
photon with respect to the final-state electron direction [[|. Note that Eq. [TJ] includes a 
factor of 1/2 to account for the identical photons in the final state. 



B. Virtual Corrections 

The matrix element for the process e~(s)j(X) — ► e~(s')j(X') is expressed by Tsai, DeR- 
aad, and Milton in the following form H, 

M%(s, sO = ^4E <0<*W ■ 4 " eA(g)«(p, (15) 
where j = 0, 1 labels the order of the matrix element and the six gauge-invariant, singularity- 



free, kinematic- zero-free, Dirac tensors £j are defined by Bardeen and Tung [ 10 1 . The authors 



calculate the six invariant matrix elements M$ within the framework of Schwinger source 
theory to order-a (M- ) and to order-a 2 (M^). They explicitly consider the case that 
the initial- and final-state electrons are longitudinally polarized and express the matrix 
element as a set of six helicity amplitudes (due to the charge-conjugation and time-reversal 
symmetries, only six of the eight matrix elements defined in Eq. [15] are independent) which 
are linear combinations of the six invariant matrix elements. The helicity amplitudes are 
then used to derive an order-a 3 expression for the unpolarized cross section which is found 



to agree with the calculation of Brown and Feynman. This cross check was found to be 
useful in locating three typographical sign errors in the rendering of the helicity amplitudes 
(which, given the complexity of the expressions, is a remarkably small number). The errors 
are listed in Appendix [A]. 

As was mentioned in the introduction, Gongora and Stuart supply spinor-product ex- 
pressions for the six tensors, 

t aa'( s > s ') = u(p r , s')ey(q') ■ & ■ e x (q)u(p, s), (16) 

which permits the application of the virtual corrections contained in the invariant matrix 
elements to the case of general initial-state and final-state electron spin directions. 

To make use of these, the system of six equations which define the helicity amplitudes was 
inverted to extract the Mj. 

The order-a 2 and order-a 3 cross sections for the process e~(s) r y(X) — > e~7 are then 
calculated (in the cm-frame) from the order-a and order-a 2 matrix elements as follows, 

^-^- ^ + 2 E ^ + m gKk^Of (IT) 

where the order-a 3 cross section has been labeled as cr^) to explicitly indicate that Eq. 
describes virtual corrections only. 



C. Soft-Photon Corrections 

The order-a 3 cross section defined in Eq. [18| contains a term that depends logarithmically 
upon a small, but nonzero, fictitious photon mass (m 7 ) used to regulate singularities in the 
virtual corrections. This unphysical term is cancelled by a similar term which arises in the 
cross section for e~7 — > e~77 for slightly massive photons. Brown and Feynman discuss this 
point at some length in Ref. 0. Since events with additional photons of energy less than 
some small value /c™ n are experimentally indistinguishable from the two-body final state, 
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they explicitly integrate the three-body cross section over the extra photon momenta q" in 
the region m 7 < E" < k™ m . The resulting soft-photon cross section is approximately equal 
to the product of a function J and the order-a 2 cross section. The order-a 3 2-body cross 
section can now be defined function of 



b 7 ' 



d 2 a (1) d 2 a ilv) r^T" d 5 a^ 

i^L (s , A; = i^ (s , A; mj) + jfj ^^(., A) 

- ~ln^ {s > A; m7) + J{m ~" k ™ m > W-w-t*' A) ' (19) 

where a^j is independent of m 7 . One should note that although Eq. |TJ is independent of 
reference frame, the actual integration was performed in the rest frame of the initial- 
state electron. The use of the resulting expression for J implies that the quantity k™ m is 
defined in that frame. 



D. The e'e+e' Final State 

The cross section for the process e~(s) r y(X) — ► e~(s')e + (s)e~(s") is calculated using 
the massive spinor-product techniques given in Ref. |J. Since the work described there 
doesn't involve positrons, its authors did not define massive positron spinors. It is extremely 
straightforward to do this from Eqs. |5] and [7] by making the replacement m — > —m. This 
interchanges the massless momentum vectors, p\ <-> p 2 , and yields the following massive 
positron spinors, 

v (p } s ) = S+ ( pl > P2 K + (p 2 ) + U _(p0 (20) 
m 

v(p, s) = - S - {P ^ P2) u + {p 2 ) + «_( Pl )- (21) 

These spinors have the correct normalization and orthogonality properties, 

u(p, s)u(p, s) = 2m, v(p, s)v(p, s) = —2m, 
u(p, s)v(p, s) = v(p, s)u(p, s) = 0. 
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Similarly, the evaluation of a massive spinor with a negative momentum leads to the re- 
placements {pi,p 2 } "~ > { — Vii —Pi} an d using Eq. [LC] one finds the correct behavior to within 
extra phases, 

v(~P, s) = iu(p, s), v(-p, s) = iu(p, s), (22) 
u{— p, s) = iv(p, s), u(—p, s) = iv(p, s). (23) 

These extra phases do not occur in the case of external photons (see Eq. [13]) and must be 
treated with some care. To avoid disturbing the phase relationships between diagrams, a 
calculation must be formulated so that all diagrams contain the same number of momentum- 
reversed massive spinors. 

The actual calculation was carried out by calculating spinor-product expressions for two 
Feynman amplitudes, D\\ and D 2 \, for the process 7 — ► e~e + e~e + as shown in Fig. |2|. These 
expressions are listed in Appendix |TJ. The matrix element is the sum of the eight diagrams 
generated by reversing one of the positron momenta and by interchanging final state electron 
momenta (according to Fermi-Dirac statistics), 

M\(s, s, s', s") = Dix(-p, s; p, s; p', s'; p", s") - D ix (-p, s; p, s; p", s"; p', s') 

- Dix{p, s; -p, s; p', s'; p", s") + D ix (p, s; -p, s; p", s"; p', s') 
+ D 2 \{-p, s;p,s;p ,s ;p ,s ) - D 2 \{-p, s;p,s;p ,s ;p,s) 

- D 2X {p, s; -p, s; p', s'; p", s") + D 2X (p, s; -p, s; p", s"; p', s'), (24) 

where: p and s are the momentum and spin of the initial-state electron, p and s are the 
momentum and spin of the final-state positron, p' and s' are the momentum and spin of 
one final-state electron, and p" and s" are the momentum and spin of the other final-state 
electron. This particular formulation does not benefit from algebraic simplifications due to a 
clever choice of the photon auxiliary momentum q. It is possible to reduce the total number 
of terms in the matrix element from 112 to 96 by defining four amplitudes instead of two 
and by choosing q appropriately. As formulated, the choice of q is truly arbitrary. 
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The matrix element given in Eq. [24] is converted into a cross section with an expression 
that is very similar to Eq. [H], 

d 5 aill , 1 



E |-Ma( S ,3,/,s")| 2 



(25) 



dE> e <m> e dE»d<l)>l v ~' ' v 64(2vr)5£ 7 (£ e + P e , _ s 
where Eg' is the energy of the second electron and (j)" is the azimuth of the second electron 
with respect to the first electron direction [0] . 



III. IMPLEMENTATION 

The Fortran-code COMRAD consists of three weighted Monte Carlo generators: 
COMTN2, COMEGG, and COMEEE. These perform integrations of the cross sections 
for the e~7, e~77, and e~e + e~ final states, respectively. Operational details are given in 
Appendix 0. 

A. Weighting Scheme 



Each of the generators produces events that consist of momentum four-vectors of the 
final-state state particles in the laboratory frame. Each event is also accompanied by a 
vector of four event weights Wj. The event weights are defined as follows: 

1 



Wi 

W 3 
W 4 



2pW(x) 
1 

1 

2pW{x) 
1 

2p(")(x) 



8,-) + -^—(s,+) 



dx n 
d n a^ 
dx n 

dx n 



dx 11 
d n a^ 
dx n 



dx n 
d n a^ 



s,+) 



(26) 
(27) 
(28) 
(29) 



dx n dx T 

where n is the dimensionality of the integrated space (n = 2 for the e~7 final state, and 
n = 5 for the three-body final states) and p- n ' is the density of trials in that space. 

The sums of the weights yield partly or fully integrated cross sections. It is convenient 
to define the following notation for these sums, 
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4°V) = E^i (so) 

<W) = X>2* (31) 
^iV)=EW3 ( 32 ) 

j 

<#W) = £W& (33) 

i 

where the variables x' define the kinematical binning chosen for a particular problem. The 
sums of the W\ and W 3 weights yield the order-a 2 and order-a 3 unpolarized cross sections, 
respectively. The sums of the W 2 and W 4 weights yield the order-a 2 and order-a 3 polarized 
cross sections and depend upon the initial spin direction s (all cross sections are given in 
millibarns). It is also convenient to define notation for the fully corrected cross sections and 
for the asymmetry functions, 

a u (x')=a£\x>)+a£\x>) (34) 

a p (s;x') = a^(s;x') + al J 1 \s;x') (35) 
0.(0) („. T f\ 

A^\s-,x') = a -^l (36) 

(37 » 

Note that the polarized cross sections are chosen to be the differences of the negative- 
helicity photon cross sections and the positive-helicity photon cross sections. In particle 
physics terminology, these are called left-handed-helicity and right-handed-helicity photons, 
respectively. One should note that in optics terminology, a negative-helicity photon is called 
Right- Circularly-Polarized (RCP) and a positive-helicity photon is called Left- Circularly- 
Polarized (LCP). 

The generation of multiple weights per event trial allows the user to significantly improve 
the statistical power of a given set of Monte Carlo trials. The uncertainty on any function 
of the four quantities Oj = {oi°\ a£\ u^} is always smaller when generated with 
correlated weights than separate, uncorrected calculations would yield. The correct estimate 
of the statistical uncertainty on any such function requires that the user accumulate the full 
4x4 error matrix Ejk, 
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ntrial 

E jk = E WjW*k, (38) 

i=l 

where the sum is over all event trials. This matrix must then be propagated correctly to the 
final result. As an example, consider the calculation of the uncertainty on the quantity, AA, 
which is the difference of the full, order-a 3 -corrected polarized asymmetry and the order-a 2 
asymmetry, 

o-(o) _|_ ^(l) o-(o) 

AA(s; x>) = A(s; x') - A® (s; x') = p {0) + ^ - (39) 

CTu ~\~ CTu CTu 

The correct uncertainty on this quantity is given by the expression, 

r/A ^ * dAA dAA . , 

S ^=^ E ^ (40 » 

where j and k label the four cross sections. 



B. Cross Checks 

The code COMRAD has been checked in a number of ways. The order-a 3 unpolarized 
cross section calculated from the unpolarized initial-state by C0MTN2 is numerically 
identical to the one calculated from the diagnostic expression given by TDM in Ref. [0 and 
to one given by Brown and Feynman in Ref. It is verified that this cross section is 
rigorously independent of the value chosen for the photon mass m 7 . It is also verified that 
the order-a 3 polarized cross section is invariant under helicity-flips of both incident particles 
(as required by parity invariance). 

The hard-photon cross section calculated by COMEGG is verified to be independent of 
the choice of photon auxiliary momenta. The polarized hard-photon cross section is found 
to be invariant under helicity-flips of both incident particles. The dependence of the cross 
sections a^j and crj^ on k™ in is shown in part (a) of Fig. |] for the case of a 50 GeV electron 
colliding with a 2.34 GeV photon (one of the cases considered by Veltman in Ref. 0]). 
Note that the cross-sections vary by approximately 1.4 mb as fc™ n is varied from 30 eV 
to 10 KeV. The sum of the cross sections, ct^\ is shown in part (b) of the figure and is 
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constant at 0.002 mb level until k™ m reaches several percent of the maximum photon energy 
and the two-body approximation for the soft-photon cross section begins to fail. Even then, 
the 10 KeV point differs by only 0.012 mb from the 30 eV point. 

The e~e + e~ cross section calculated by COMEEE is found to be independent of the 
choice of photon auxiliary momentum. The polarized e~e + e~ cross section is found to be 
invariant under helicity-flips of both incident particles. 

The sum of the virtual, soft-photon, and hard-photon cross sections calculated by COM- 
RAD is compared with the numerical result presented in Ref. Q for the case of a 50 GeV 
electron colliding with a 2.34 GeV photon. The ratio of the unpolarized cross sections 
cr^/a^ is presented as a function of the laboratory energy of the scattered electron E' lab 
in part (a) of Fig. |j. The COMRAD calculation predicts that a^' / increases from from 
—0.14% near the kinematical edge at 17.90 GeV to +0.2% near the beam energy. The Velt- 
man calculation predicts that the ratio decreases from +0.3% near the edge to +0.2% near 
the beam energy. The physically correct behavior follows from a simple kinematical analysis. 
In the center-of-mass frame, the emission of an additional photon reduces the energy and 
momentum available to the scattered electron. Given the large mass of the electron, the 
fractional change in the momentum P' e is larger than than fractional change in the energy 
E' e . The laboratory energy of a backscattered electron is given by the following expression, 

E[ ab = 7 (E' e - P' e ) , (41) 

where 7 is the Lorentz factor for the highly-boosted cm-frame (the velocity is assumed to be 
one). It is straightforward to show that although E' e and P' e are decreased by the emission 
of an additional photon, the difference E' e — P' e increases. The laboratory energy of the 
backscattered electron is therefore increased by the emission of an additional photon. It is 
clear that photon emission depopulates the Compton kinematical edge region and that 
should be negative near the endpoint. 

The order-a correction to the longitudinal polarization asymmetry is shown in part (b) 
of Fig. |. The quantity AA(s z ; E[ ah ) (defined in Eq. |39|) predicted by COMRAD is compared 
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with the similar quantity given in Ref. ||]. Good agreement is observed. 



IV. RESULTS 

This section describes the application of the COMRAD code to several accelerator- 
based polarimetry cases. The first case (the SLD polarimeter) deals with the detection of 
the scattered electrons to measure longitudinal polarization. The second case (the HERA 
polarimeters) involves the detection of scattered photons to measure longitudinal and trans- 
verse electron polarization. The final Linear Collider polarimeter) illustrates the 
detection of final state electrons when there is sufficient energy to produce the e~e + e~ final 
state. 



A. The SLD Polarimeter 

The SLD Polarimeter |J is located 33 m downstream of the SLC interaction point (IP). 
After the 45.65 GeV longitudinally-polarized electron beam passes through the IP and before 
it is deflected by dipole magnets, it collides with a 2.33 eV circularly-polarized photon beam 
produced by a pulsed frequency-doubled Nd:YAG laser. The scattered and unscattered 
components of the electron beam are separated by a dipole-quadrupole spectrometer. The 
scattered electrons are dispersed horizontally and exit the vacuum system through a thin 
window. A multichannel Cherenkov detector observes the scattered electrons in the interval 
from 17 to 27 GeV/c. 

The helicities of the electron and photon beams are changed on each beam pulse ac- 
cording to pseudo-random sequences. Each channel of the Cherenkov detector measures 
the asymmetry in the signals Sj observed when the electron and photon spins are parallel 
(| J z \ = 3/2) and anti-parallel (| J z \ = 1/2), 

i ~ ^(3/2) + ^-(1/2) ~ e 7 J K ] 
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where: j labels the channels of the detector, VI is the electron beam polarization, P 7 is the 
photon polarization, and Aj is the analyzing power of the j channel. The analyzing powers 
are defined in terms of the Compton scattering cross section and the response function of 
each channel Rj, 

where E[ ab is the laboratory energy of the scattered electron. The unpolarized cross sec- 
tion a u (E[ ah ) and longitudinal polarization asymmetry A(s z ; E[ ah ) are shown as functions of 
scattered electron energy in Fig. ^. The order-a 2 quantities are shown as dashed curves in 
parts (a) and (b) of the figure. The fractional correction to the unpolarized cross section 
a^/a^ is shown as the solid curve in part (a) of the figure. It increases from —0.2% near 
the endpoint at 17.36 GeV to +0.2% at the beam energy. The correction to the asymmetry 
function AA(s z ; E[ a A is shown as the solid curve in part (b) of the figure. Near the endpoint, 
AA is almost exactly one thousand times smaller than the order-a 2 asymmetry. It becomes 
fractionally larger near the zero of A<® at 25.15 GeV. 

The effects of the order-a 3 corrections upon the analyzing powers of the seven active 
channels of the Cherenkov detector are listed in Table [II The nominal acceptance in scattered 
energy, the order-a 2 analyzing power Af\ and the order-a 3 fractional correction to the 
analyzing power are listed for each channel. The SLC beam polarization is determined 
from the channels near the endpoint (5-7). It is clear that proper inclusion of the radiative 
corrections increases the analyzing powers by 0.1% of themselves. This decreases the beam 
polarization by the same fractional amount. Since the left-right asymmetry is the ratio of 
the measured Z-event asymmetry Az and the beam polarization, 

' e 

the application of the order-a 3 corrections increases the measured value of Alr by 0.1% 
of itself. The corrections are much too small and have the wrong sign to account for the 
SLD/LEP discrepancy. 
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B. The HERA Polarimeters 



The e ± ring of the HERA e^-p collider is the first storage ring to operate routinely 
with polarized beams and it is the first storage ring to operate with a longitudinally polarized 
beam [O] . The ring is instrumented with transverse and longitudinal Compton polarimeters. 



1. The Transverse Polarimeter 



The HERA transverse polarimeter [12| collides 2.41 eV photons from a continuous-wave 
Argon-Ion laser with the 27.5 GeV HERA positron beam. The scattered photons are sep- 
arated from the electron beam by the dipole magnets of the accelerator lattice and are 
detected by a segmented tungsten-scintillator calorimeter located about 65 m from the e + -7 
collision point. The scattering rate is sufficiently small that the calorimeter measures the 
energy and vertical position of individual photons. 

When the positron beam is transversely polarized, the differential cross section depends 
upon the azimuthal directions of the scattered particles. The average direction the scattered 
photons changes when the laser helicity is reversed. The polarimeter measures the projected 
vertical direction 9 y and energy k' lah of each scattered photon. The shift in the centroid 
of the 9 y distribution that occurs with helicity reversal S9 y neas (k[ ah ) is proportional to the 



product of the photon polarization and the vertical positron polarization V\ 



y 

e ' 



^r aS ( fc U) = (0 V )- - (0y} + = VI ■ V, ■ S9 y (k[ ah ), (45) 

where S9 y (k[ ah ) is the shift for 100% positron and photon polarizations. This quantity is 
given by the following expression, 

se y (kU) = g/^vW,)^™^ (46) 

°"«(^lab) 

where 9' and <f)'^ are the polar angle and azimuth of the scattered photon in the laboratory 
frame. Note that 9' is a constant for fixed k[ ah . 
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The order-a 3 corrections modify the function 59 y {k[ ab ). The exact modification depends 
upon the details of how the polarimeter reacts to the two-photon final state. It is assumed 
that the segmented calorimeter of the HERA transverse polarimeter cannot distinguish 
between one-photon and two-photon final states. The energy measured by the calorimeter 
for two-photon final states is then the sum of the individual photon energies k' lah = k[ ah (l) + 
k{ ah (2). The measured vertical angle is assumed to the the energy- weighted mean of the 
individual photon angles, 9 y = [0 y (l)A;{ ab (l) + y (2)k[ ah (2)] /k[ ah . The order-a 2 function 
59 y °\k' lah ) is plotted as function of laboratory photon energy in Fig. |6|. The maximum 
angular separation of 5.6 jum occurs near 8 GeV. The fractional change caused by the order- 
a 3 corrections A59 y /59 y °* > is also shown as a function of k' lah . Note that the correction is 
typically +0.08% near the maximum separation which would lower the measured transverse 
polarization by the same fractional amount. 



2. The Longitudinal Polarimeter 



A longitudinal polarimeter at HERA has been built by the HERMES Collaboration |]13 
The 27.5 GeV HERA positron beam is brought into collision with a 2.33 eV photon beam 
produced by a pulsed frequency- doubled Nd:YAG laser. The scattered photons are separated 
from the electron beam by the dipole magnets of the accelerator lattice and are detected by 
an array of NaBi crystals. Since several thousand scattered photons are produced on each 
pulse, it is not possible to measure the cross section asymmetry as a function of photon 
energy. Instead, the calorimeter measures the asymmetry in deposited energy as the 
photon helicity is reversed, 

rpdep rpdep 

Ae = Edep ~ E t p = V:V,A E , (47) 

where E± p is the energy deposited by all accepted photons in the crystal calorimeter. The 
analyzing power Ae is given by the following expression, 

/ dk{ ah k[ ah R(k[ ah )a u (k[ ah ) 



where R(k[ ah ) describes the response of the detector. For this estimate, it is assumed that the 
calorimeter has uniform response in energy from the minimum accepted energy of 56 MeV 
(lower energy photons miss the calorimeter) to the maximum energy of 13.62 GeV. The 
order-a 2 analyzing power and the full order-a 3 correction are 

A%> = 0.1838 (49) 
A - 4 (0) 

E (0 ) E = +0.20%. (50) 
A E 

The fractional correction to the longitudinal polarization scale is therefore —0.20%. 



C. A Linear Collider Polarimeter 

Longitudinally polarized beams are likely to be important features of a future Linear 
Collider. It is assumed that any such machine will include SLC-like polarimetry which 
detects and momentum-analyzes scattered electrons. The unpolarized cross section a u (E[ ah ) 
and longitudinal polarization asymmetry A(s z ; E[ ah ) are shown as functions of scattered 
electron energy in Fig. [7] for the case of a 500 GeV electron beam colliding with a 2.33 eV 
photon beam. The order-a 2 quantities are shown as dashed curves in parts (a) and (b) 
of the figure. The cross section is largest near the backscattering edge at 26.42 GeV. The 
longitudinal asymmetry function is 0.9944 at the kinematic endpoint. It decreases rapidly 
with increasing energy and passes through zero at 50.19 GeV. The fractional correction to 
the unpolarized cross section a^/a^ is shown as the solid curve in part (a) of the figure. 
It increases from —1.6% near the endpoint to +1.2% at the beam energy. Superimposed 
upon this is the contribution of the e~e + e~ final state which is kinematically constrained to 
the region 34.36 GeV < E[ ah < 386.1 GeV. The effect of this final state is to increase the 
correction to the 1.0-1.7% level in the kinematically allowed region. The correction to the 
asymmetry function AA(s z ; E^A is shown as the solid curve in part (b) of the figure. Near 
the endpoint, AA is — 4xl0 -4 and represents a negligible correction. Due to the influence 
of the e~e + e~ final state, it decreases to — 2.2xl0 -3 near 49 GeV then begins to increase to 
+5.3xl0 -3 near 306 GeV where it is a 1% correction to the asymmetry function. 
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V. SUMMARY 



The construction of a computer code, COMRAD, to calculate the cross sections for the 
spin-polarized processes e~7 — > e~7, e~77, e~e + e~ to order-a 3 has been described. The 
code is based upon the work of Tsai, DeRaad, and Milton [0 for the virtual and soft- 
photon corrections. The hard-photon photon corrections and the application of the virtual 
corrections to arbitrary electron spin direction are based upon the work of Gongora and 
Stuart [PJ. The calculation of the cross section for the e~e + e~ final state was performed 
by the author. As implemented, the code calculates cross sections for circularly-polarized 
initial-state photons and arbitrarily polarized initial-state electrons. Final-state polarization 
information is not presented to a user of the code but is present at the matrix element level. 
The modification of the code to extract this information would not be difficult. 

The order-a 3 corrections to the longitudinal polarization asymmetry calculated by COM- 
RAD agree well with those of Veltman ||]. However, the order-a 3 corrections to the unpo- 
larized cross section calculated by COMRAD do not agree with those of Veltman. 

The application of the code to the SLD Compton polarimeter indicates that the order- 
a 3 corrections produce a fractional shift in the SLC polarization scale of —0.1%. This 
shift is much too small and of the wrong sign to account for the discrepancy in the Z-pole 
asymmetries measured by the SLD Collaboration and the LEP Collaborations. 

The application of the code to the photon-based polarimeters at the HERA storage 
ring indicates that the order-a 3 corrections also have small effect on the measurements of 
the HERA positron polarization. The effects on the transverse polarization measurements 
are typically less than 0.1%. The effect upon the calibration of the HERMES longitudinal 
polarimeter is a somewhat larger 0.2%. 

The application of the code to a polarimeter at a future Linear Collider indicates that 
the order-a 3 corrections are very small near the Compton edge but increase to the 1% level 
elsewhere. The e~e + e~ final state contributes significantly to the net corrections. 
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APPENDIX A: ERRATA 

The following typographical errors were found in Ref. ||: 

1. All of the spinor products of the form u±(qi)$u T (q 2 ) given in Eqs. 3.3-3.10 and C.1-C.8 
formally vanish (they are the traces of an odd number of gamma matrices) and should 
be replaced by u±(qi)p'u±(q 2 ) (the helicities of the w-spinors are correct in all cases 
and the helicities of the w-spinors are wrong in all cases). The right-hand-sides of the 
spinor product definitions are nearly all correct (see items 4 and 5 below). 

2. The sign of the second term in square brackets on the right-hand-side of Eq. 3.5.1 in 
the preprint is incorrect [— s + (p 2 , q") • • • should be +s + (p 2 , q") • • •]. Unfortunately, a 
serious typesetting error in the published version significantly altered the equation. 
The correct equation should read as follows, 

s+(pi,P2)s-(px,q) 



D ++ -{q,P2]q',P2]q",Pi) 



x 



s+(P2, q)s+(p 2 , q')s-(pi, q") 

s+{p'i,q")s-{p'i,p'2) - , w / n- / /w / \ 
5 MPi¥fe M -G°2 )u- [q )p a u-{p 2 ) 

+ s + (P2 , q" ) s - (pi , p' 2 ) u _ ( q' )i> a u _ (p 5 



3. The sign of the fourth term in square brackets on the right-hand-side of Eq. 3.6.1 is 
incorrect [— s + (p 2 , q') ■ ■ ■ should be +s + (p 2 , q') ■ ■ ■}. 

4. The right-hand-side of the first of Eqs. 3.6.2 is incorrect. The quantities p\ and p 2 
should be replaced by p[ and p' 2 , respectively. 
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5. The left-hand-side of the second of Eqs. 3.7.2 is incorrect. The quantity p[ should be 
replaced by p 2 . The right-hand-side is also incorrect. The sign of the second term 
should be flipped [— S-(q",p' 2 ) ■ ■ ■ should be +S-(q" ,p' 2 ) ■ ■ ■]. 

6. The second factor in the second term in square brackets on the right-hand-side of 
Eq. 3.8.1 should be S-(pi,q") instead of S-(p[,q"). 

7. The fourth factor in the fourth term in square brackets on the right-hand-side of 
Eq. C.7.1 should be S-(p' 2 ,q') instead of s_(p' 2 ,q). 

8. The heading of Eqs. 4.10 which states that they define quantities of the form e'_£je + 
is correct and all of the left-hand-sides which state the reverse helicity configuration 
are wrong. 

9. The heading of Eqs. 4.11 which states that they define quantities of the form e+£je_ 
is correct and all of the left-hand-sides which state the reverse helicity configuration 
are wrong. 

The following typographical errors were found in Ref. 0: 

1. The signs of two of the tree-level helicity amplitudes given in Eqs. 5 are incorrect. The 
signs of the amplitudes — h; H — ) [the third amplitude] and f^ 2 \ — h; ++) [the fifth 
amplitude] should be reversed. 

2. The sign of the order-a 2 amplitude — h; H — ) defined in Eq. 8 should also be 
reversed. 

APPENDIX B: THE MATRIX ELEMENT FOR £T 7 -» E'E^E' 

The matrix for the process e~(s)j(X) — > e~(s')e + (s)e~(s") is calculated from the two 
amplitudes for the process 7(A) — * e + (s)e~~ (s')e + (s)e~ (s") shown in Fig. |2|. This formulation 
is chosen so that each term in Eq. ^4] has exactly one negative momentum. The internal 
momenta shown in the Fig. |2| are defined as, 
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Q = p" + P 

Pa = ~ + Q) = V - Q 

Pb = q- p = P+ Q- 
Using the techniques described in Ref. and the Chisholm identity 

7 /1 M±(gi)7/^±(g 2 ) = 2 [u±(q 2 )u±(qi) + u^{q 1 )u^{q 2 )\ , 



(Bl) 
(B2) 
(B3) 



(B4) 



it is straightforward to evaluate the amplitudes D\\ and D 2 \. Unfortunately, the exact form 
for each of these functions depends upon the initial-state photon helicity A. They are listed 
below: 

2\/2e 3 f 



Di+(p,s;p,s;p',s';p",s") = — , 

-2p' ■ q Q 2 s + {q,q) { 

S-(p' 2 ,q)s+(pi,p 2 )s + (q,p 1 )s-(p 2 ,p 2 ) + s + (g,p , 1 )s_(p' 1 ,p , 2 )s_(g,p , 2 / )s + (pi,pi) 

+S-(p / 2 ^) s +(Pi ) Pi)[ s +(^Pi) s -(p'i ) P2) + s+(q,p 2 )s-(p' 2 ,p 2 ) - s+(q,q)s-(q,p 2 ') 
s+(q,p'i)s-(p' 1 ,p' 2 )s + (p 1 ,p 2 )s-(p , 2 \p 2 ) l _ 



MK>P 2 ) g +(P2,Pl) 

S-(p 2 ^)s+(pi,p 2 )s + (g,Pi)s_(p2,P2) + s + (g,p , 1 )s_(p' 1 ,p , 2 )s_(g,p 2 )s+(p 1 ,pi) 

+s-(p , 2^) s +(Pi ) Pi)[ s +(^K) s -(Pi ) P2) + s+(<?,P 2 )s-(i4P2) - s+(g,g)s_(g,p 2 ) 
, s + (q,p' 1 )s-(p' 1 ,p , 2 )s + (pi,p 2 )s-(p2,P2) 



S-(q,Pi)s+{p[,p") + s-{q,p' 2 )s + {p' 2 ,p") > (B5) 



2v^e 3 



Di-(p,s;p,s;i/,8 , ;]/ , ,8") = — , 

-2p' • g Q z s-(q, q) [ 

+s-(p 2 ,g)s+(pi,Pi)[s + (g ) pi)s-(pi,P2 / ) + p' 2 )s-(j4p 2 ) 
, s+(g,pi)s_(pi,j^)a+(pi,p2)s-(p?,P2) r 



s-(g,Pi)s+(pi,Pi) + s_(g,j£)s+(p£,pi) 



+ 



s-(pi,p 2 )s+{p 2 ,p 1 ) 



s -(P 2 i q)s+(pi,p 2 )s+(q, p")s_(p 2 , p 2 ) + s + (g,p / 1 )s_(p' 1 ,p , 2 )s_(g,p 2 )s + (p' 1 / ,pi) 
+S-(P2,?) s +(Pi,Pi) s +(97P'i) s -(p'i>P2) + s+(g,p 2 )s_(p 2 ,p 2 ) 
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+S-(q,P2)s+(p' 2 ,Pi) - s-(q,q)s+(q,p") 



S-(q,p[)s + (p[ : Pi) 



2y/2e 3 



D2+(p,s;p,s;p ,s ;p ,s ) , , ■ 

-2p- g Q 2 s+(g,g) [ 

-s+ (g, pi ) s_ (p' 2 , P2 ) [s+ (pi , pi) s- (pi , g) + s+ (pi , p 2 ) s_ (p 2 , g) 

s-(p 2? 9)s+(9,g) - s-{p 2 ,pi)s+(pi,q) 



-S-(j%,P2)s+(p2,q) 



+ 



S~(P'1,P2) S +(P2,P1) 



-s_(p / 1 ,p 2 )s + (g,Pi)s + (p' 1 ,Pi) s _(p 2 ,g) + s + (pi,p 2 )s_(g,p 2 )s_(p / 2 ,p 2 )s + (p / 1 / ,g) 

-S-(P2,p2)s+(q,pi)[s + (p",p 1 )s-(p 1 ,q) + s+(p",p 2 )s-(p 2 ,q) 
s+(p' 1 ,Pi)s^(p' 1 ,p' 2 )s + (p u p 2 )s_(q,p 2 ) 



s_(p 2 ,g)s+(g,g) 



-s_(p 2 ,Pi)s+(pi,g) - s_(p 2 ,p 2 )s+(p 2 ,g) 



2>/2e 3 



A>_ (p, s;p,s; p', s'\ p", s") = — . ^ , 

-2p • g Q 2 s_(q, q) { 

-S-(pi,P 2 )s+(?,Pi)s+(p'i,Pi)s-(p 2 ,?) + s+(Pi,P2)s-(g,P 2 )s-(p 2 ,P 2 ) s +(Pi,9) 

+s + (g,Pi)s_(p 2 ,p 2 )[s + (pi,g)s_(g,g) - s+(p~i,pi)s_(pi, g) - s+(pi,p 2 )s_(p 2 , g) 

, s+(p / 1 ,pi)s-(pi,p 2 )s + (pi,p 2 )s_(g,p 2 ) r / // x / v, / // \ / v 
+ 2 s_ (p 2 ,Pi)s+(pi,q) + s-{p 2 ,p 2 )s + {p 2 ,q) 

S -(Pl,P2) S +(P2,Pl) f 



-s_(p / 1 ,p 2 )s+(g,Pi)s + (p' 1 ,Pi) s _(p 2 ,g) + s + (pi,p 2 )s_(g,p 2 )s_(p 2 ,p 2 )s + (pi,g) 

+s-(p 2 ,P2)s+(g,Pi)[s+(Pi,g)s-(g,g) - s+(K,pi)s-(pi,g) - s+(p",p 2 )s_(p 2 ,g) 
. s+(pi,K) s -(K>P2) s +(pi>P2)s-(g,p 2 ) r v / x, v / v 

+ 2 s_(p2,Pi)s+(pi,g) + s_(p 2 ,p 2 )s+ p 2 ,g) 
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APPENDIX C: OPERATIONAL DETAILS 



The Fortran-code COMRAD [T3| consists of a main program COMRAD that controls 
the three weighted Monte Carlo generators: COMTN2, COMEGG, and COMEEE. These 
perform integrations of the cross sections for the e~7, e~77, and e~e + e~ final states, re- 
spectively. All three generators use a common set of conventions, input parameters, and 
a common interface routine called WGTHST. The routine WGTHST permits the user to 
accumulate event weights in a manner that is appropriate to his/her needs. Note that all 
quantities discussed in this section are assumed to be of type REAL*8 unless otherwise 
specified. 



1. The Program COMRAD 

The program COMRAD initializes all quantities and sequentially calls each of the event 
generators. Communication with the generators occurs through the / CONTROL / common 
block which is specified within COMRAD. This common block contains the variables: EB, 
EPHOT, XME, XMG, KGMIN, ALPHA, PI, ROOT2, BARN, SPIN(3), LDIAG, LBF, and 
NTRY. 

The variables EB and EPHOT specify the energy and laboratory frame to be used in the 
calculation. It is assumed that the incident electron is moving in the +z-direction with an 
energy EB GeV (EB > m). The incident photon is assumed to be moving in the — z-direction 
with an energy EPHOT GeV. The spin of the initial-state electron is specified in its rest 
frame by the three- vector SPIN(3). The maximum energy of additional soft-photons and the 
minimum energy of hard-photons is also specified in the electron rest frame {k™ in ) by the 
variable KGMIN. The integer variable NTRY sets the number of trials for each of the event 
generators (COMTN2 and COMEGG generate NTRY trials whereas COMEEE produces 
smaller event weights and generates NTRY/20 trials). The logical flags LDIAG and LBF 
activate the calculation of the Tsai-DeRaad-Milton and Brown-Feynman expressions for the 
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corrections to the unpolarized cross section in COMTN2 (for diagnostic purposes). The 
common block /CONTROL/ also contains several constants used by the generators. 



2. The Generator COMTN2 

The subroutine COMTN2 simulates the two-body e~7 final states. The calculation is 
carried out in the center-of-mass frame using Eq. [19] and Eqs. |26| - p9| to calculate the event 



weights VF1-W4. The density function p^ 2 \Q!^) is chosen to be uniform in the polar variables 



cos6*g and 4>' e , 



P (2) m = ^f, (CI) 



where N tria i is the number of event trials. 



3. The Generator COMEGG 



The subroutine COMEGG simulates the three-body e 77 final states. The calculation 



is carried out in the center-of-mass frame using Eq. [14] and Eqs. ^ and |29] to calculate the 
event weights W 3 and W4 (W\ and W 2 are always returned as 0). The five quantities E' e , 
cos9' e , 0g, E' , and 7 are chosen according to the density function as follows, 



Ntri 



trial 



An ■ 2n E™ ax + E min — E' e 



7 7 



I fTmin CV 



(C2) 



where: N tria i is the number of generated trials (some are later discarded); E™ ax and E™ ax are 
the maximum electron and photon energies in the cm-frame; E™ in is the minimum photon 
energy in the cm-frame; and C e and C 7 are normalization constants given as follows, 

1 



a 



111 



max 1 T?min 

7 



m)/E, 



mm 

7 



2 In 



E™ ax /E, 



mm 

7 



(C3) 
(C4) 



The minimum energy in the cm-frame is related to the minimum photon energy in the 
initial-state electron rest frame as follows, 
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Ef n = — — k™ in , (C5) 

-l- J cm 

where E cm is the total center-of-mass energy. Note that a photon emitted in the cm-frame in 
the —z direction with energy E™ m has an energy /c™ m in the electron rest frame. If emitted 
in any other direction, it has a smaller energy in the electron rest frame. 

After all five variables have been chosen, the electron and photon energies are checked 
for consistency with three-body kinematics (the angle 9 ei between the electron and photon 
directions must satisfy the condition |cos# e7 | < 1). If this condition is not satisfied, the 
trial is discarded. If it is satisfied, the four-vectors p', q', and q" are generated. The photon 
energies in the initial-state electron rest frame are then calculated and if either is found to be 
less than k™ m , the trial is discarded. The kinematical boundary of the integration is therefore 
exactly the same as the one that defines the upper-limit of the soft-photon integration (and 
the function J). The integrated region is thus the complement of the soft-photon region and 
the sum of the cross sections returned by COMTN2 and COMEEG is independent of k™ in . 
The event generation procedure retains approximately 60% of the generated trials over a 
wide range incident electron and photon energies. 



4. The Generator COMEEE 



The subroutine COMEEE simulates the three-body e e + e final states. The calculation 



is carried out in the center-of-mass frame using Eq. g5| and Eqs. ^ and |29] to calculate the 
event weights W 3 and W4 {W\ and W 2 are always returned as 0). The five quantities E' e , 
cos9' e , 0' e , E", and <\>" e are chosen according to the density function p^> as follows, 

where N tr iai is the number of generated trials (some are later discarded) and E™ ax is the 
maximum electron energy in the cm-frame. Note that the use of uniform phase space works 
well near threshold (the polarimetry case) but is inadequate at very high energies. 
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After all five variables have been chosen, the electron energies are checked for consistency 
with three-body kinematics (the angle e / e « between the electron directions must satisfy the 
condition |cos0 e / e //| < 1). If this condition is not satisfied, the trial is discarded. This 
procedure retains approximately 70% of the generated trials near the e~e + e~ threshold. 



5. The Interface Routine WGTHST 

The routine WGTHST allows the user to accumulate the information needed for his/her 
purposes. The routine is called by the main program once before any event generation to 
permit initialization. It is called by each of the event generators (COMTN2, COMEGG, 
and COMEEE) at the end of each event trial. And finally, it is called by the main program 
after return from the last generator to permit the information to be output. 

All communication with the routine occurs through the argument list, 

SUBROUTINE WGTHST(IFLAG,NEM,PP,NGAM,QP,NEP,PB,WGT), 

where: IFLAG is an integer flag which indicates the initialization call (0), an accumulation 
call (1), or the output call (2); NEM is an integer which indicates the number of electrons 
in the final state (1 or 2); PP(4,2) contains the four-vectors of the NEM electrons in the 
laboratory frame; NGAM is an integer which indicates the number of photons in the final 
state (0-2); QP(4,2) contains the four- vectors of the NGAM photons in the laboratory 
frame; NEP is an integer which indicates the number of positrons in the final state (0 or 
1); PB(4) is the four- vector of the positron; and WGT(4) contains the four weights W\-W± 
defined in Eqs. |26|-|29]. Note that the event weights have been defined such that correctly 



normalized total cross sections are obtained by summing the weights exactly once per call to 
WGTHST. The calculation of final-state particle yields therefore requires that the weights 
be accumulated each time the given type of particle is encountered. 
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DEC-OpenVMS operating systems. 
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FIGURES 
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FIG. 1. Amplitude for the process e — > e 777 (time flows left to right). 
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FIG. 2. Two amplitudes for the process 7 — > e e + e e + (time flows left to right). 
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50 GeV e~on 2.34 eVy 
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FIG. 3. The unpolarized virtual plus soft-photon cross section a£j and hard-photon cross 
section fTg'yy Oil"© shown in part (a) as functions of /c™ n for the case of a 50 GeV electron colliding 
with a 2.34 eV photon. The sum of these cross sections 0$ is shown in part (b). 
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50 GeV e~ on 2.34 eV y 
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FIG. 4. The ratio Ou /<7„ is shown in part (a) as a function of the laboratory energy of the 
scattered electron for the case of a 50 GeV electron colliding with a 2.34 eV photon. The difference 
in the fully order-a 3 corrected asymmetry and the order-a 2 asymmetry is shown as a function of 
the laboratory energy of the scattered electron in part (b). The COMRAD calculation is shown as 
solid lines and the calculation of Ref. [Ml is shown as dashed lines. 
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45.65 GeV e - on 2.33 eVy 
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FIG. 5. The order-a 2 cross section and the ratio /<Ju^ are shown in part (a) as functions 
of the laboratory energy of the scattered electron for the case of a 45.65 GeV electron colliding 
with a 2.33 eV photon. The order-a 2 asymmetry function and the difference in the fully order-a 3 
corrected asymmetry and the order-a 2 asymmetry are shown as functions of the laboratory energy 
of the scattered electron in part (b). 
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FIG. 6. The order-a 2 function S9y'\k[ ab ) is plotted as the dashed curve as a function of the 
laboratory energy of the scattered photon for the case of a 27.5 GeV positron colliding with a 
2.41 eV photon. The fractional difference in the fully order-a 3 corrected function and the order-a 2 
function is shown the solid curve. 
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500 GeV e - on 2.33 eV y 
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FIG. 7. The order-a 2 cross section and the ratio / are shown in part (a) as functions 
of the laboratory energy of the scattered electron for the case of a 500 GeV electron colliding with a 
2.33 eV photon. The order-a 2 asymmetry function and the difference in the fully order-a 3 corrected 
asymmetry and the order-a 2 asymmetry are shown as functions of the laboratory energy of the 
scattered electron in part (b). The vertical dashed-dotted lines indicate the allowed kinematic 
region for electrons from the e~7 — ► e~e + e~ subprocess. 
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TABLES 

TABLE I. The effect of order— a 3 radiative corrections upon the analyzing powers of the SLD 
Compton polarimeter. 



Channel £ e Acceptance (A-A^)/A^(%) 

7 17.14-18.02 GeV 0.7133 0.096 

6 18.02-19.00 GeV 0.6483 0.097 

5 19.00-20.11 GeV 0.5520 0.103 

4 20.11-21.38 GeV 0.4309 0.118 

3 21.38-22.83 GeV 0.2851 0.153 

2 22.83-24.53 GeV 0.1228 0.285 

1 24.53-26.51 GeV -.0396 -.673 
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